Viscosity in Strongly Interacting Quantum 
Field Theories from Black Hole Physics 



P.K. Kovtim, 1 * D.T. Son, 2 t and A. 0. Starinets 3 * 

Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 
2 Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA 
3 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 

December 2004 



Abstract 

The ratio of shear viscosity to volume density of entropy can be used to characterize 
how close a given fluid is to being perfect. Using string theory methods, we show that 
this ratio is equal to a universal value of h/AirkB for a large class of strongly interacting 
quantum field theories whose dual description involves black holes in anti-de Sitter 
space. We provide evidence that this value may serve as a lower bound for a wide class 
of systems, thus suggesting that black hole horizons are dual to the most ideal fluids. 
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Introduction. — It has been known since the discovery of Hawking radiation pQ that black 
holes are endowed with thermodynamic properties such as entropy and temperature, as first 
suggested by Bekenstein [2] based on the analogy between black hole physics and equilib- 
rium thermodynamics. In higher-dimensional gravity theories there exist solutions called 
black branes, which are black holes with translationally invariant horizons For these 
solutions, thermodynamics can be extended to hydrodynamics — the theory that describes 
long- wavelength deviations from thermal equilibrium jl]. In addition to thermodynamic 
properties such as temperature and entropy, black branes possess hydrodynamic character- 
istics of continuous fluids: viscosity, diffusion constants, etc. From the perspective of the 
holographic principle [SHE], a black brane corresponds to a certain finite-temperature quan- 
tum field theory in fewer number of spacetime dimensions, and the hydrodynamic behavior 
of a black-brane horizon is identified with the hydrodynamic behavior of the dual theory. 
For these field theories, in this Letter we show that the ratio of the shear viscosity to the 
volume density of entropy has a universal value 

- = — %- « 6.08 x 1(T 13 Ks. (1) 
s 4irk B 

Furthermore, we shall argue that this is the lowest bound on the ratio rj/s for a wide class 
of thermal quantum field theories. 

Viscosity and graviton absorption. — Consider a thermal field theory whose dual holo- 
graphic description involves a .D-dimensional black-brane metric of the form 

ds 2 = g [ $ N dx M dx N = /(0 (dx 2 + dy 2 ) + g^) d^dC • (2) 

[The 0(2) symmetry of the background is required for the existence of the shear hydrody- 
namic mode in the dual theory, thus making the notion of shear viscosity meaningful.] One 
can have in mind, as an example, the near-extremal D3-brane in type IIB supergravity, dual 
to finite-temperature M = 4 supersymmetric SU(iV c ) Yang-Mills theory in the limit of large 
N c and large 't Hooft coupling 00131101, 
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but our discussion will be quite general. All black branes have an event horizon [r = ro for 
the metric ©], which is extended along several spatial dimensions [x, y, z in the case of 
The dual field theory is at a finite temperature, equal to the Hawking temperature of the 
black brane. 

The entropy of the dual field theory is equal to the entropy of the black brane, which is 
proportional to the of its event horizon, 

s =4 w 
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where G is Newton's constant (we set h = c = kg = 1). For black branes A contains a trivial 
infinite factor V equal to the spatial volume along directions parallel to the horizon. The 
entropy density s is equal to a/ (AG), where a = A/V. 

The shear viscosity of the dual theory can be computed from gravity in a number of 
equivalent approaches JTJ^lEI]. Here we use Kubo's formula, which relates viscosity to 
equilibrium correlation functions. In a rotationally invariant field theory, 

v= l im _L / dbdxe** ({T xy (t,x), 2^(0,0)]) . (5) 

Here T xy is the xy component of the stress-energy tensor (one can replace T xy by any com- 
ponent of the traceless part of the stress tensor). We shall now relate the right-hand side 
of (0) to the absorption cross section of low-energy gravitons. 

According to the gauge-gravity duality [TU] , the stress-energy tensor T^ u couples to metric 
perturbations at the boundary. Following Klebanov let us consider a graviton of 

frequency uj, polarized in the xy direction, and propagating perpendicularly to the brane. In 
the field theory picture, the absorption cross section of the graviton by the brane measures 
the imaginary part of the retarded Greens function of the operator coupled to h xy , i.e., T xy , 

a abs H = -^!lmG» = - f dtdxe™ ([T xy (t,x), T xy (0,0)]) , (6) 
u uj J 

where k = V8nG appears due to the normalization of the graviton's action. Comparing 
and (jBJ), we find 

r] = ^^ = l^G- (7) 

Graviton absoprtion cross section at low energy. — The absorption cross section <7 a b s is 
calculable classically by solving the linearised wave equation for h y . We now show that h x y = 
h xy / f obeys the equation for a minimally coupled massless scalar in the background (J2J). This 
is similar to cosmological tensor perturbations on a Friedmann-Robertson- Walker spacetime, 
which obey the equation for a massless scalar field [15] . 

Consider small perturbations around the metric, Qmn — 9mn + ^mat- We assume that 
the only non- vanishing component of Hmn is h xy , and that it does not depend on x and 
V'- h xy = h xy (£). This field has spin 2 under the 0(2) rotational symmetry in the xy 
plane, which implies that all other components of Iimn can be consistently set to zero |12j . 
Einstein's equations can be written in the form 

T 

Rmn = Tmn — — — - Qmn > (8) 

where the stress-energy tensor Tmn and its trace T depend on other fields such as the 
dilaton and various forms supporting the background (J2J), for example, the fields appearing 
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in the low-energy type II string theory Again, 0(2) xy rotational symmetry implies that all 
perturbations of matter fields can be set to zero consistently. Thus when M and N are x or 
y, the right-hand side of Einstein's equations reads (a, j3 — x, y) 

T ( T(°) \ 

T aP - - — - g a( 3 = -[C + - — - (8 af3 f + hap) , (9) 



D-2 a ^ V D-2 t 

where C is the Lagrange density of matter fields and is the trace of the unperturbed 
stress-energy tensor. 

T M n = 9m n £ H , (10) 

where £ represents the Lagrangian density of the fields, and dots denote terms of second 
and higher orders in the perturbation h xy . Substituting the unperturbed metric © into 
Einstein's equations, one finds 



□/ (df) 



T(°) 

C+ —2- <"> 



/ P 

Expanding Einstein's equations to linear order in h xy , one has 

1 1 (df) 2 ( T(°) \ 

R xy = —-Dh xy + -&>f - K -^h xy = -{£+ — h xy . (12) 



Combining Eqs. and (|12p . we obtain an equation for h 



xy 



(Of) 2 nf 

nh xy -2-±- d^h xy + 2 h xy j- h xy — . (13) 



Changing the variable from h xy to = h xy /f, one can see that h y indeed satisfies the 
equation for a minimally coupled massless scalar: Oh y = 0. The absorption cross section of 
a graviton is therefore the same as that of the scalar. 

The absorption cross section for the scalar is constrained by a theorem ^3 QB| > which 
states that in the low-frequency limit uj — > this cross section is equal to the area of the 
horizon, <7 a b s = gl Since s = a/AG, one immediately finds that 

2 -i4-- (14 » 

s Attkb 

where h and ks are now restored. This shows that the ratio rj/s does not depend on the 
concrete form of the metric within the assumptions of Refs. f7[ QB] • Indeed, this ratio is the 
same for Dp- ([TTJ El), M2- and M5- (|19j) branes and for deformations of the D3 metric 
^21l2ni- This fact is very surprising, given that the corresponding dual field theories are very 
different. We do not have an explanation for the constancy of rj/s in these theories based on 
field-theoretical arguments alone. 
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A viscosity bound conjecture. — Most quantum field theories do not have simple gravity 
duals. Is our result relevant in a broader setting? We speculate that the ratio rj/s has a 
lower bound 



for all relativistic quantum field theories at finite temperature and zero chemical potential. 
The inequality is saturated by theories with gravity duals. 

One argument supporting the bound (JT5j) comes from the Heisenberg uncertainty prin- 
ciple. The viscosity of a plasma is proportional to er m f t , where e is the energy density and 
r m f t is the typical mean free time of a quasiparticle. The entropy density, on the other hand, 
is proportional to the density of quasiparticles, s ~ A^n. Therefore rj/s ~ kg r^e/n. Now 
e/n is the average energy per particle. According to the uncertainty principle, the product 
of the energy of a quasiparticle e/n and its mean free time r m f t cannot be smaller than h, 
otherwise the quasiparticle concept does not make sense. Therefore we obtain, from the 
uncertainty principle alone, that n/s > h/ks, which is (115)) without the numerical coefficient 
of l/(47r). We also conclude that rj/s is much larger than h/ks in weakly coupled theories 
(where the mean free time is large). 

Another piece of evidence supporting the bound (JT5j) comes from a recent calculation [21] 
of rj/s in the J\f = 4 supersymmetric SU(N C ) Yang-Mills theories in the regime of infinite 
N c and large, but finite, 't Hooft coupling g 2 N c . The first correction in inverse powers of 
g 2 N c corresponds to the first string theory correction to Einstein's gravity. The result reads 



where £(3) ~ 1.2020569... is Apery's constant. The correction is positive, in accordance with 
(115)1 . It is natural to assume that r)/s is larger than the bound for all values of the 't Hooft 
coupling (Fig. 1). 

The bound ([15)1 . in contrast to the entropy bound [22] and Bekenstein's bound [23J, 
does not involve the speed of light c and hence is nontrivial when applied to nonrelativistic 
systems. However, the range of applicability of (|15[) to nonrelativistic systems is less certain. 
On the one hand, by subdividing the molecules of a gas to an ever-increasing number of non- 
identical species one can increase the entropy density (by adding the Gibbs mixing entropy) 
without substantially affecting the viscosity. On the other hand, the conjectured bound is far 
below the ratio of rj/s in any laboratory liquid. For water under normal conditions rj/s is 380 
times larger than h/(4nkB)- Using standard tables [2UI2S] one can find rj/s for many liquids 
and gases at different temperatures and pressures. Figure 2 shows temperature dependence 
of r]/s, normalized by K/ (kirks), for a few substances at different pressures. It is clear that 
the viscosity bound is well satisfied for these substances. Liquid helium reaches the smallest 




(16) 
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Figure 1: The dependence of the ratio rj/s on the 't Hooft coupling g 2 N c in M = 4 su- 
per symmetric Yang-Mills theory. The ratio diverges in the limit g 2 N c — > and approaches 
h/A-irks from above as g 2 N c — > oo. The ratio is unknown in the regime of intermediate 
't Hooft coupling. 

value of r)/s, but this value still exceeds the bound by a factor of about 9. We speculate that 
the bound (fT3|) is valid at least for a single-component non-relativisitic gas of particles with 
spin or 1/2. 

Discussion. — It is important to avoid some common misconceptions which at first sight 
seem to invalidate the viscosity bound. Somewhat counterintuitively, a near-ideal gas has a 
very large viscosity due to the large mean free path. Likewise, superfluids have finite and 
measurable shear viscosity associated with the normal component, according to Landau's 
two-component theory. 

The bound (|15jl is most useful for strongly interacting systems where reliable theoretical 
estimates of the viscosity are not available. One of such systems is the quark-gluon plasma 
(QGP) created in heavy ion collisions which behaves in many respects as a droplet of a 
liquid. There are experimental hints that the viscosity of the QGP at temperatures achieved 
by the Relativistic Heavy Ion Collider is surprisingly small, possibly close to saturating the 
viscosity bound [26J. Another possible application of the viscosity bound is trapped atomic 
gases. By using the Feshbach resonance, strongly interacting Fermi gases of atoms have 
been created recently. These gases have been observed to behave hydrodynamically 27 and 
should have finite shear viscosity at nonzero temperature. It would be very interesting to test 
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Figure 2: The viscosity-entropy ratio for some common substances: helium, nitrogen and 
water. The ratio is always substantially larger than its value in theories with gravity duals, 
represented by the horizontal line marked "viscosity bound." 

experimentally whether the shear viscosity of these gases satisfies the conjectured bound. 

This work was supported by DOE grant DE-FG02-00ER41132, the National Science 
Foundation and the Alfred P. Sloan Foundation. 
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